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I. INTRODUCTION 

Extensions of the scalar sector of the Standard Model (SM) beyond the minimal single Higgs 
doublet are of great interest in model building and collider phenomenology and are, as yet, largely 
unconstrained by experiment. Such extensions arc common in models that address the hierarchy 
problem of the SM, such as supersymmetric models [1] and little Higgs models [21], as well as in 
models for neutrino masses, dark matter, etc. Most of these extensions contain additional SU(2)l- 
singlet, -doublet, and/or -triplet scalar fields. However, some extensions of the SM contain scalars 
in larger multiplets of SU(2)l. Such larger multiplets have been used to produce a natural dark 
matter candidate [3HS] , which is kept stable thanks to an accidental global symmetry present in the 
Higgs potential for multiplets with isospin T > 2. Three different models with a Higgs quadruplet 
(T = 3/2) have also been proposed for neutrino mass generation 6 8 . Models in which the SM 
SU(2)-doublet Higgs mixes with a seven-plet (T = 3), aided by additional representations of SU(2), 
have also been studied recently in Ref. [§] ■ 

In this paper we consider models that extend the SM scalar sector through the addition of a single 
large multiplet. For multiplets with n = 2T + 1 > 5 (isospin 2 and larger), the scalar potential 
of these models always preserves an accidental global U(l) or Zi symmetry at the renormalizable 
level. If unbroken, such a symmetry forces the lightest member of the large multiplet to be stable. 
Spontaneous breaking of an accidental global U(l) symmetry is phenomenologically unacceptable 
because it would lead to a massless Goldstone boson that couples to fermions through its mixing 
with the SM Higgs doublet's neutral Goldstone, and thus mediate new long-range forces between 
SM fermions. Furthermore, perturbative unitarity of scattering amplitudes involving pairs of scalars 
and pairs of SU(2) gauge bosons requires that T < 7/2 (i.e., n < 8) for a complex scalar multiplet 
and T < 4 (i.e., n < 9) for a real scalar multiplet [TTj] . 

The models that preserve such an accidental global symmetry can be grouped into three classes 
based on the hypercharge Y of the large multiplet, as follows: 

(i) Models with a real, Y = multiplet, with n = 5, 7, or 9, corresponding to isospin 2, 3, or 4 
(a real multiplet must have integer isospin). The large multiplet is odd under an accidental 
global Z2 symmetry. These models have been considered in Refs. [3l 0] as possible candidates 
for "next-to-minimal" dark matter. 

(ii) Models with a complex multiplet with n = 5, 6, 7, or 8, with Y = 2T (we work in the 
convention Q — T 3 + Y/2). The large multiplet is charged under an accidental global U(l) 
symmetry. The hypercharge is chosen so that the lightest member of the multiplet can be 
electrically neutralj^j The masses of the states in the large multiplet are split by an operator 
of the form (&T a &)(X^T a X), where $ is the SM Higgs doublet, X is the large multiplet, 
and r a and T a are the appropriate SU(2) generators. We study these models in the current 
paper. 

(Hi) Models with a complex multiplet with n = 6 or 8, with Y = 1. The large multiplet is odd 
under an accidental global Z2 symmetry. The would-be accidental global U(l) symmetry 
is broken by an operator of the form ($^r°$)(X^T a X), where <&, X denote the conjugate 
multiplets. Such an operator can appear only when n is even. We will address these models 
in a forthcoming paper |12) . 



Models in which the lightest member of the large multiplet is electrically charged are excluded or strongly con- 
strained by the absence of electrically-charged relics. Metastable multi-charged states are constrained by direct 
collider searches to be heavier than about 400-500 GeV, depending on their charge 
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In this paper we study the constraints on the models with a complex multiplet with n = 5, 6, 
7, or 8, with Y = 2T. We first determine the constraints on the scalar quartic couplings from 
perturbative unitarity and precision electroweak measurements. We then examine the bounds on 
the neutral scalar \° from cosmological considerations. Our goal is not to determine whether x° is 
able to account for the entire observed quantity of dark matter this possibility is strongly excluded 
by dark matter direct detection experiments — but rather to evaluate the ultimate viability of the 
model as a target for collider searches. Assuming a standard thermal history of the universe, we 
compute the thermal relic density of x° an d X°* ancl compare it with the limits from dark matter 
direct detection experiments. In conjunction with the requirement that m x o > mz/2 ~ 45 GeV 
from the invisible width of the Z boson, we find that these cosmological constraints exclude the 
models with n = 6, 7, and 8. 

The n — 5 (T — 2) multiplet has a dimension-5 Planck-suppressed interaction with the SM Higgs 
doublet of the form 

C D — ^$$$$A f + hue, (1) 
Mpi 

where Mpi is the Planck mass. This operator induces a mixing of the SM Higgs into the neutral 
component of X, \' = X°' r ~ ^4>°' r , with e ~ v3 /[{ m x (> ~ m \) M Pi]- Here v ~ 246 GeV is the 
SM Higgs vacuum expectation value. The neutral member of X can then decay via its Higgs 
component. Assuming that \° is the lightest state, its lifetime ranges from a few days to a few years 
for m x o ~ 100-1000 GeV. This puts the decays of the lightest neutral state of the n = 5 model well 
after big-bang nucleosynthesis and well before the recombination surface of the cosmic microwave 
background radiation. Direct detection constraints from present-day experiments therefore do not 
apply to this model, and it remains viable. 

The paper is organized as follows. In Sec. [Il] we give the Lagrangians and mass spectra for the four 
models that we consider. In Sec.|III|we obtain the indirect constraints on the model parameters from 



unitarity and the oblique parameters, and comment on collider constraints. In Sec. IV we calculate 
the upper bound on the relic density of neutral scalars x ) X°* from dark matter direct-detection 
experiments. In Sec. [V] we compute the relic density from thermal freeze-out and show that in 
all cases it yields a relic density too large to be consistent with the bound from direct detection. 



We conclude in Sec. VI Feynman rules, formulas for the oblique parameters, and formulas for the 



partial decay width of the SM Higgs to two photons are collected in the appendices. 



II. THE MODELS 



The models that we consider extend the SM through the addition of a single complex scalar 
multiplet X, with hypercharge Y — 2T (normalized so that Q = T 3 +Y/2). The hypercharge is 
chosen so that the lightest member of X can be neutral. The size of the multiplet X is restricted 
to n = 2T + 1 < 8 by the requirement that tree-level amplitudes for SU(2) gauge bosons scattering 
into the states in X, W a W a — > x®*X® \ remain perturbative [TU]. When n > 5, the scalar potential 
possesses an accidental global U(l) symmetry corresponding to phase rotations of X. This U(l) 
symmetry ensures that the lightest member of X is stable, at least at the level of renormalizable 
operators. These two conditions leave us with four models to consider, with n = 5, 6, 7, and 8. 

The gauge-invariant scalar potential is given by, 



V($, X) = m 2 $ f $ + M 2 X^X + Ai ($ f $) 2 + \ 2 &$XiX + A 3 $ t r a $V t T a V + C(A 4 ), (2) 
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where r a and T a are the generators of SU(2) l in the doublet and n-plet representations, respectively. 
$ is the usual SM Higgs doublet, and the large scalar multiplet takes the form 

x = (x +M V-,x ) T 



The mass of particle \® with charge Q = T 3 + Y/2 > is given by, 



m 2 Q 



M 2 + -v 2 



\o — — X'lT 
i 2 



M 2 + -v 2 



A 2 - -A 3 ( Q 



(3) 



(4) 



where v ~ 246 GeV is the SM Higgs vacuum expectation value (vev). The neutral particle x°, with 
T 3 = -T = -{n - l)/2, will have a mass, 



m 2 x o = M 2 



A 2 + ~A 3 (n-l) 



= M 2 + -v 2 A , 



(5) 



where we define the dimensionless coupling Ao as the quantity in brackets above. The masses of 
the charged states x can be written in terms of the x° mass as, 



1 



m x Q = m x o 



v 2 X 3 C 



(6) 



We require that the stable lightest member of X is electrically neutral; this forces us to take A3 < 0. 



III. CONSTRAINTS ON COUPLINGS AND MASSES 

A. Unitarity constraints on scalar quartic couplings 

The scalar quartic couplings A 2 and A3 given in Eq. ([2]) can be bounded by requiring perturbative 
unitarity of the zeroth partial wave amplitude. The partial wave amplitudes are related to scattering 
matrix elements according to 

M = 167r^(2J+l)a. / F J (cos6»), (7) 
J 

where J is the orbital angular momentum of the final state and Pj(cos9) is the corresponding 
Legendre polynomial. Perturbative unitarity of the zeroth partial wave amplitude dictates the 
tree- level constraint, 

|Rea | < \. (8) 

The coupling A 2 controls the strength of the isospin-zero x*X ~^ amplitude, while A3 con- 
trols the strength of the isospin-one x*X ~ * 4>* $ channel. Because we are working with large scalar 
multiplets, the isospin-zero x*X ~ * WW,BB and isospin-one x*X ~ * WB amplitudes can be sig- 
nificant [Til], leading to more stringent coupled-channel limits on A 2 and A3. We neglect all other 
contributing processed] and work in the high-energy limit. 



2 Additional contributions to the matrix of coupled-channel amplitudes come from quartic couplings of X as well 
as (p*(j> <f>*<f> amplitudes proportional to Ai. We find numerically that including these contributions generically 
leads to a slightly tighter constraint on A2 and A3, but this constraint depends on the interplay between the Ai 
contributions and those from the quartic X couplings. 
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The relevant amplitudes for the isospin-zero channels are, 



a (lx*x]o -> [M]o) = -g^^, 



a ([x*x]o [W] ) 
a ([x*x]o -► [BB]o) 



g 2 (n 2 — l)y/n 
16tt 271 ' 

2 167rc?„ 2 



16?r c 2 ^ 



(9) 



where the \*X ~> WW,BB amplitudes include both of the contributing transverse gauge boson 
polarization combinations [10] and we used Y = 2T = n — 1 in the last line. Here g is the SU(2)l 
gauge coupling and sw, cw = sin 6yy, cos 0w are the sine and cosine of the weak mixing angle. We 
define the following normalized isospin-zero field combinations, 



'Jo = 



[x*x]o 
[WW]o 



1 n— 1 
. /n Z — / 



X Q *X Q , 



2=0 



1 

75 

[BB] = BB/V2 



V2W + W~ 



(W 3 W 3 



(10) 



The relevant amplitudes for the isospin-one channels are, 
a o([x*x]i -> 



yjn{n 2 



1) 



32^ 



A, 



ao([x*x]i -> = 



Y^n(n 2 - 1) ff 2 s w (n - l)^/™^ 2 - 1) 



167T Cw 



76 



(11) 



167r -v/6 

where again the x*X ~ * WB amplitude includes both of the contributing transverse gauge boson 
polarization combinations |10j . Here we used the following normalized isospin-one field combina- 
tions, 



1 



T2 y 



[x*xh 



12 



n(n 2 



12 



2=o 



n(n 2 



1) ? 



Q* 



2=0 



Q 



n — 1 



[WB]i = W 3 B. 



(12) 



Finding the largest eigenvalue of each coupled-channel matrix and applying the unitarity con- 
straint of Eq. (18]), we find the unitarity bounds on A2 and A3, 



|A 2 | < 



/32tt 2 
n 



An 2 -!) 2 
1 24 



5 4 T 
c w 



(n-iy 



|A,| < 2, 



/ 384tt 2 
n(n 2 — 1) 



s 2 
r 2 



{n-lf 



(13) 
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n 


lAiri 


|a!ti 


5 


7.64 


n.i 


6 


6.49 


8.17 


7 


5.01 


6.11 


8 


2.17 


4.41 



TABLE I. Upper limits on | A2 1 and | A3 1 from perturbative unitarity, for Y = 2T 



Recall that A3 must be negative so that x° is the lightest member of the l arge multiplet. A2 can 
have either sign. Numerical bounds are given for n = 5, 6, 7, and 8 in Table W 



B. Electroweak precision constraints 

The multiplet X contributes to electroweak observables through the oblique parameters S, T, and 
U [T3]. The contributions of such a scalar multiplet obeying a U(l) global symmetry, so that the 
mass eigenstates have definite T 3 , were computed for arbitrary isospin and hypercharge in Ref. |14j . 
We summarize the results in Appendix [B] for completeness. 

The contributions to the oblique parameters depend only on n, Y, and the masses ui x q of each 
state in the multiplet. Therefore, for a given n and setting Y — 2T = n — 1, the oblique parameters 
constrain only two model parameters, which can be chosen as m x o (which sets the overall mass 
scale) and A3 (which controls the mass splittings). 

The current experimental values relative to the SM with Higgs mass = 126 GeV are S exp = 
0.03 ± 0.10, T cxp = 0.05 ± 0.12, f7 cxp = 0.03 ± 0.10, with relative correlations of p ST = 0.89, 
Ptu — —0.83, psu — —0.54 [15] • We use these values to constrain m x o and A3 via a two-parameter 
X 2 variable; for details see Appendix [b| 95% confidence level (% 2 = 5.99) limits are shown for the 
two parameters m x o and Am = rn x + — m x a in Fig. [TJ 

At low m x o the constraint is dominated by the S parameter and leads to an upper bound on Am 
that is linear in m x o . This bound can be parameterized as 



Am 



in 



,+ — m x o = 0.031 



n — 4 



0.13 



(14) 



For larger m x o ~ 1 TeV, the constraint from the T parameter becomes important and limits the 



value of Am independent of m x o. For m^ 



5-6 TeV, the unitarity limit on A3 becomes the 



strongest constraint on Am, as shown by the dashed lines in Fig. [T] 

Notice that for n > 6 the mass splitting between x° and the next-lightest state x + 1S constrained 
to be no more than about 12 GeV; for m x o ~ 100 GeV the splitting is less than 1.5 GeV. The 
maximum allowed mass splitting decreases with increasing n. 



In this table we use g 



,2 
, W 



■■ 0.231, and a = 1/128. 




FIG. 1. 95% confidence level constraints on Am = m x +i — m x o as a function of m x o from the S, T, and 
U parameters, for the scalar multiplets with n = 5, 6, 7, and 8 and Y = 2T = n—l. Dashed lines indicate 
the upper limit on Am from the unitarity bound on A3. The left panel shows the low-m x o region while the 
right panel extends to higher masses. 



C. Direct collider constraints 



Scalar particle masses below about 100 GeV are constrained by x®X®* P a i r production in e + e~ 
collisions at the CERN Large Electron-Positron (LEP-II) collider. However, a dedicated search for 
the decay signatures in the models we consider here has not been made; furthermore, x^X®* events 
may be difficult to detect if the mass splittings are small, leading to low-energy charged particles 
from the x® decays. 

Regardless of these difficulties, the LEP-I measurement of the Z boson invisible decay width 
puts a stringent constraint on Z — > X°X°* independent of the mass splittings. This leads to the 
requirement 

m x o > mz/2 ~ 45 GeV. (15) 

As we will show, this constraint together with cosmological considerations is sufficient to exclude 
the models with n = 6, 7, and 8. Further investigation of the direct collider constraints on the 
n = 5 model is beyond the scope of this paper. 



IV. DARK MATTER DIRECT DETECTION CONSTRAINT 



The allowable relic density of X° is constrained by its non-observation in direct dark matter 
detection experiments. In particular, the fraction of the ambient dark matter density that can be 
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attributed to \ is bounded from above according to 



O -limit 

X <^V, (16) 



where (Jgj is the spin-independent, per-nucleon scattering cross section for \° 01 X°* an d ^gf^ * s 
the experimental upper limit on the spin-independent dark matter scattering cross section obtained 
assuming the canonical ambient dark matter density. The strongest experimental upper limit 
currently comes from the XENON100 experiment [IB]. We note that this limit will not apply to 
the model with n — 5 because in this model the relic \° particles decay away on a timescale short 
compared to the age of the universe. 

Because x° is a complex scalar, it scatters off nucleons via both Z and Higgs exchange. The Z 
exchange diagram yields a large scattering cross section, leading to very stringent constraints on 
fl x /fluM- We compute the cross sections in the zero- velocity limit assuming equal densities of x° 
and such that f2 x = O x o + fi x o». We have, 

where af nt is the interference between the Z and Higgs exchange diagrams. For x we have 

0. 0. 

a z h = a z h an d a int — ~ a iat> s0 tne interference term cancels in the total scattering rate for equal 
densities of x° an d X°* • 

The Z-exchange cross section for scattering off a single nucleon N = p, n is given by 



where v ~ 246 GeV is the SM Higgs vev, m^r is the nucleon mass, and the vector couplings of the 
Z to the nucleon are given by the sum of the corresponding valence quark couplings|j 



v 2m z ( 1 2 \ j, v 2m z ( 1 



The axial- vector couplings do not contribute in the zero- velocity limit. Notice that a\ becomes 
independent of m x o in the large-m x o limit, where it is of order m 2 N /v A . Note also that er^ is fixed 
with no free parameters once m x o and the size of the multiplet n are specified. 

The ^.-exchange cross section for scattering off a single nucleon N — p, n is given by 

where A v is the hx°X°* coupling defined in Eq. ^ and the Higgs-nucleon Yukawa couplings are 
given by [15] 



fn = — (0.350 ± 0.048), = — (0.353 ± 0.049). (21) 

v v 



4 These nucleon vector couplings do not receive any QCD corrections in the limit of zero momentum transfer due 
to the conservation of the vector current |17| . 
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Notice that goes like l/m^ in the large-m x o limit, where it is of order m%/v m^ . The Higgs- 
exchange contribution is thus generically much smaller than the Z-exchange contribution. The 
Higgs-exchange contribution also depends on the parameter A = A 2 + Aa(n — l)/4. We obtain the 
least stringent upper bound on n x /fluM when Ao = 0. 

Because dark matter particles moving in the galactic halo have de Broglie wavelengths that are 
large compared to the size of a nucleus, the amplitudes for scatt erin g off e ach nucleon add coherently. 



This can be accounted for by replacing f^- and /jy in Eqs. (18) and (20 1 above by the coherent 
nucleon-averaged values, 



{f V N? {f V N f 



[Zfl + (A-Z)fi] 2 
A 2 



u h N ? -+ u h N f = [zfhp + {A A2 z)f!l]2 - u';? - if h n )\ (22) 

where Z is the atomic number and A the atomic mass of the nucleus. For xenon, Z = 54 and 
A ranges from 124 to 136. We make a weighted average over the natural abundances of xenon 
isotopes [T§] . 

The upper limit on S} x /£Idm from XENON100 [16] as a function of the x° mass is shown in 
Figs. [2] and [3] for the multiplets with n = 6, 7, and 8 (the limit depends on n like l/(n — l) 2 )- The 
shaded area bounded by the red curve is excluded. Here we have set Ao = in order to obtain the 
most conservative limit; taking Ao ^ has only a tiny effect on the limit. 

Finally we comment on the behavior of the upper limit on f2 x /SlDM at large \° masses. The 
scattering cross section <7gj is overwhelmingly dominated by the Z exchange contribution, which is 
independent of m x o in the large-m x o limit. The XENON100 collaboration quotes a cross section 
limit for dark matter masses up to 1000 GeV. When the dark matter particle mass is large compared 
to the mass of the target nucleus, the energy transfer for a given target nucleus asymptotes to a 
constant which depends only on the velocity of the incoming dark matter particle. The experimental 
cross section limit then varies inversely with the ambient number density of dark matter particles, 
which in turn goes like the (fixed) mass density times l/m x o. Therefore, the upper bound on 
f2 x /fiDM grows linearly with m x o for masses large compared to the target nucleus mass, and can 
be extrapolated to arbitrarily heavy masses. 



V. THERMAL RELIC DENSITY 



We now compute the thermal relic density of x° + X°* and compare it to the upper bound 
from direct detection found in the previous section. We assume a standard thermal history of the 
universe: i.e., that the temperature was high enough at one time for x° to have been in thermal 
equilibrium, that there were no late-decaying relics that would significantly dilute the abundance 
of x ^*''; a nd that there were no late-decaying relics that decayed to states in X and hence boosted 
the x ^*'' relic abundance. 

We will show that, for all allowed parameter choices, the thermal relic abundance of x ^ is too 
large to be consistent with the upper bound on 57 x /J7dm from direct detection. The U(l)-preserving 
models with n = 6, 7, and 8 are thus excluded assuming a standard thermal history. 

The fraction of the dark matter density that is due to x is given in terms of the total x°X°* 
annihilation cross section by 

_ (ciy e ;)std ^ 23 ^ 

°dm \(<?v re i(x°X a * -> any)) ' 
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jq-8 I , , ^Aj , , ^_^J , , , I , , l_^J , , , I , . . 

10 50 100 500 1000 5000 lxlO 4 

m x , (GeV) 



FIG. 2. The fraction f2 x /f2oM of the total dark matter density for the n = 6 model as a function of m x o. 
The shaded area above the red curve is excluded by direct detection constraints from XENON100 data |16| . 
conservatively taking Ao = 0. The solid blue curve shows the predicted relic density assuming a standard 
thermal history of the universe, for the parameters Ao = 0.01, A3 = —0.01. The black dashed curve is the 
relic density in the case that co-annihilations are maximal (see text for details). Masses below rriz/2 (to 
the left of the vertical black line) are excluded by the LEP constraint on the invisible Z width. 




10 50 100 500 ]000 5000 ] x if) 4 10 50 100 500 1000 5000 | x I0 4 



m x , (GeV) m x , (GeV) 

FIG. 3. The same as Fig. [2] but for n = 7 (left) and 8 (right). 

where v re i is the relative velocity of the two colliding dark matter particles and (av re i) s td is the 
"standard" annihilation cross section required to obtain the correct dark matter relic abundance, 
for which we use (cv r ei)std = 3 x 10 -26 cm 3 /s 20J. The brackets indicate an average over the 
velocity distribution at the time of freeze-out, which is only necessary if the annihilation cross 
section vanishes in the v re i — > limit. The factor of 1/2 in the denominator accounts for the 
probability that, in a collision, any given \ particle meets one with the opposite U(l) charge so 
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that an annihilation can take place. This ratio is shown by the solid blue lines in Figs. [2] and [3j 
We explain the ingredients in what follows. 



A. Annihilation to two-body final states 



A x°X°* P an " can annihilate to the two-body final states W + W~, ZZ, hh, and //. We compute 
the annihilation cross sections in the zero- velocity limit. Because these cross sections are all nonzero 
in this limit, we do not need to average over the velocity distribution. 

The annihilation cross sections to two-body final states are given in the v re i — > limit by, 






( m 2 
3-4—^- 


m x° 


V m w 


o 

777 n 

3 * 


4 




m w 


m w 


) 


' A 2 


( m l° 
3-4^ 


2 

mfn 

X 


V m z 


2 

x" 






m| 


m% j 


+ 



av re i 



m h 

2 

x° 



4m; 





9 \ ^/^ 2 A 2 

{Am 2 x o ~m 2 h ) 2 ' 




Here the coefficients in the cross sections to WW and ZZ are given by, 

v 2 A Q n 4(n-l) 



A 



w 



(n-l) + 



4m 2 o 



-4; 



(n-l) 2 + 



v 2 Ao 



4rn 2 



Bw — 



Bo 



2 2 2 

vf x o x + 



4(rt- l) 2 
2rn 2 — ml- 



(24) 



(25) 



Diagrams involving s-channel Z-exchange are zero in the low-energy limit, so that annihilation to 
// proceeds only through an s-channel Higgs. We checked our analytic results by implementing 
the relevant couplings into CalcHEP [21] . 

Above threshold, x°*X° ~^ ZZ has the largest cross section, followed by x°*X° ~^ W + W~ , which 
is smaller by about a factor of 20 for n = 6. Annihilation rates to hh and ti are much smaller, 
as can be seen by the fact that their kinematic thresholds are not even visible in Figs. [2] and [3j 
The latter two processes are controlled by the coupling Ao; we took the sample value Ao = 0.01 in 
Figs. [2] and [3j As we will see in Sec. |V C| significantly larger values of A are constrained by the 
measured rate for h — > 77. Ao also contributes to the annihilations to WW and ZZ through the 
Aw and Az coefficients; its effect is numerically small and falls with increasing m x o. We also took 
the sample value A3 = —0.01; A3 has a tiny effect on the annihilation cross section to WW through 
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the x°~X + mass splitting. Overall, the total annihilation cross section above the WW threshold 
depends very weakly on Ao and A3, and is instead controlled almost entirely by n and m x o. 

The cross sections for annihilation to WW, ZZ, and hh fall like l/m? at large m x o, while the 

cross section to // falls like l/m* . This leads to the growth of £I x /Qdm proportional to for 

the upper bound on Q x /Qdm from direct 
detection grows only linearly with m x o. Increasing m x o thus leads only to more severe conflict 
between the relic abundance and the direct detection limit. 

We finally comment on the possibility that attractive electroweak interactions between x° and 
X°* form bound states that increase the annihilation cross section for m x o 3> mw, an effect known 
as Sommerfeld enhancement [52]. This effect was studied in Ref. [23] for a real scalar 7-plet with 
hypercharge zero, which found an increase in the mass at which the dark matter candidate obtained 
the correct relic density by about a factor of 3 (to about 25 TeV) , corresponding to about an order 
of magnitude enhancement of the annihilation cross section during freeze-out. However, for our 
models, the direct detection constraint is 10 5 -10 6 times stronger than the perturbative relic density 
at large m x o > 1 TeV. We expect that an enhancement of the annihilation cross section sufficiently 
large to affect our exclusion would be extremely difficult to obtain. 



large m x o shown in Figs. As we saw in Sec. |IV| 



B. Annihilation to off-shell WW below threshold 



Below the WW threshold, the largest two-body annihilation process is x°*X° ~^ bb, which has 
a small cross section. Annihilation to off-shell WW can be significantly larger. We compute the 
annihilation cross section to off-shell WW by generating x°*X° e+v & e ~v e using CalcHEP [2T] 
and then multiplying by [1/BR(M / — > ev)] 2 = 81 at tree level. We include bb and off-shell WW in 
the blue solid curves in Figs. [2] and [3J for m x o < mw- We again take Aq = 0.01 and A3 = —0.01. 
Annihilations to off-shell ZZcovdd also be included; however, their contribution is very small 
compared to off-shell WW because the Z bosons are further off shell for a given m x o . 



C. Resonant annihilation through the Higgs pole 

The most interesting feature in the annihilation cross section below the WW threshold is the 
Higgs pole at m x o = rrih/2. The possibility that this resonant annihilation may suppress the X relic 
density enough to evade the direct detection constraints is excluded by a combination of constraints 
from the oblique parameters and the observed rate for pp — > h — > 77 from the CERN Large Hadron 
Collider (LHC) 

The x°X°* h bb annihilation cross section near the Higgs resonance is obtained from the 



last line of Eq. ( 24 ) by making the replacement in the denominator, 



(4m x o - mlf -> (4m^ - + m^. (26) 



Here is the total decay width of the SM Higgs boson; the contribution from bb final states is 
given at tree level by [2"5] . 

N c mlm h ( 4mg x 3 J 



'h 



where N c = 3 is the number of colors. 
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62.5(1 

j (GeV) 



FIG. 4. Fractional relic density Q. x /Qdm for m x o ~ mj,/2 from x°X°* — ► ^ annihilation, for various values 
of Ao (solid, dashed, and dotted black curves). Also shown are the upper limits on Q x /Qy)m from direct 
detection for (top to bottom) n = 6, 7, and 8 (red horizontal lines). 



Note that we can capture the effects of higher-order corrections and additional final states on 
c v rei{x a X°* h an y) near the Higgs resonance simply by choosing the value of mj, to yield 
the correct total SM Higgs width, T h = 4.07 ± 0.16 MeV [M] for m h = 125 GeV. We obtain this 
width from the tree-level formula above when mi, — 4.08 GeV. Using this value of mj, we compute 
S1 x /(]dm m the vicinity of the Higgs resonance using the tree-level x°X°* ~^ bb cross section and 
h ^ bb partial width. The result is shown in Fig. [4] for various values of AoQ 

The cross section for x°X°* ~^ bb is proportional to Aq, which in turn depends on A2 and A3. 
These couplings also control the one-loop contribution of x® to the Higgs decay to two photons. 
Setting to y o = m/j/2, the S, T and U parameters put an upper bound on |A 3 |, as summarized in 
Table [nj 
within t 



(A3 < is required in order for x to be the lightest member of the multiplet). For A 3 
le allowed range, we compute the partial width T(h — > 77) as a function of A2, leading to 
a constraint on Ao. The contributions of the charged x states to T(h — > 77) are summarized in 
Appendix [C] This partial width normalized to its SM value is shown as a function of A in Fig. [5] 
for the two extreme cases, A3 = and the limiting value allowed by the S, T, and U parameter 
constraints. 

We see that increasing |Ao| quickly drives T(h — > 77) to unacceptably large values. In particular, 
\A-o\ ^ 0.09 - required to evade the direct detection limit for the n — 6 model - is ruled out by 
the LHC measurement of the signal strength for pp —> h — >• 770 Allowing a very generous factor 
of three enhancement over the SM prediction (already excluded at more than 2a by both ATLAS 



5 When specified in terms of Ao, o"tV e ; (x°X°* ~~ ► bb) is independent of n. 

6 In the models considered here, the Higgs couplings are identical to those in the SM except for the contributions 
to the loop-induced /177 and h'yZ couplings from loops of charged scalars. Because the decays to 77 and 7Z 
contribute only a tiny fraction of the Higgs total width, the LHC Higgs signal strength in the 77 channel is given 
to a very good approximation by fi = cr(pp — > h — > ~n)/asM(pP — > h — > 77) — r(ft — > "ff) /T s m (h 77)- 
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n 






5 


-1.4 x 10" 


2 


6 


-6.0 x 10" 


3 


7 


-3.3 x 10" 


3 


8 


-2.0 x 10" 


3 



TABLE II. 95% confidence level lower limit on A3 from the S, T, and U parameter constraints as a function 
of the size of the multiplet, for m x o = m;,/2 = 62.5 GeV. 




FIG. 5. The partial width for h — > 77 normalized to its SM value as a function of Ao, for m x o = rrih/2 — 
62.5 GeV and A3 = (left) and the limiting (minimum) value allowed by the S, T, and U parameter 
constraints as given in Table |H| (right) . 



and CMS [24]) requires |Ao| < 0.042 for n — 6; the limit becomes more stringent for larger n. 
After applying these constraints, we see that resonant annihilation through the Higgs pole cannot 
suppress Q x /Qdm enough to allow us to evade the direct detection limits. 



D. Co-annihilations 



We finally consider the possibility that the mass splittings among the states \® are very small, 
as is favored by the oblique parameter constraints for low x° masses. In this case, all members of 
the multiplet can be present in the thermal bath during freeze-out. This opens the possibility of 
annihilations involving the electrically charged states to 77, Zj, and W ± j final statesj^] which 
are on shell below the WW threshold. Such electroweak-strength annihilation cross sections to 



7 The cross section for x®X®* ~* h"f vanishes in the zero-velocity limit. 
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two-body final states can easily dominate over those to three-body final states and the bottom- 
Yukawa-suppressed annihilation cross section to bb. 

We evaluate the potential impact of such co-annihilations by considering the extreme case in 
which all states X are present with equal abundances during freeze-out. Two conditions are 
required: the mass splittings must be very small so that the equilibrium thermal abundances of 
each species are equal, and the quartic couplings multiplying operators 0{X i ) must be large enough 
to maintain equal abundances of states with each charge, even though the annihilation cross section 
to SM states is different for states with different charge. For larger mass splittings, co-annihilations 
become less important, and the situation relaxes to our original analysis. 

In the limit of equal abundances of all states X^*\ the fraction of the dark matter density due 
to x i s given in terms of the annihilation cross sections by 

_ (ctVei)std , 2g N 

"dm l^Y, Ql ,Q 2 >o(<™rel(x Ql X Q ** ^*ny))' 

Note in particular the new factor 1/n 2 in the denominator, which represents the average over initial 
charge states for a multiplet of size n. The sum runs over all initial charge combinations; only 
the combinations with Qi — Q2 = and ±1 contribute below the WW threshold. We include 
annihilations to the two-body final states 77, Zj, W ±, y, and bb. Cross section formulas are given 
in Appendix [D] The cross sections for the first three processes are independent of A and depend 
on A 3 only through the masses m x Q. For the scalar couplings we set A = 0.01 as usual but we 
now take A 3 = 0, corresponding to degenerate masses for all 

The resulting relic abundance of \ m the full co-annihilation case is shown by the black dashed 
lines in Figs. [2] and [3] The relic abundance is dramatically reduced compared to that obtained 
considering only x°X annihilations. However, co-annihilations still do not allow us to evade the 
direct detection limits: the relic density of x remains about four orders of magnitude above the 
direct detection exclusion bound for mz/2 < m x o < mw- 



VI. CONCLUSIONS 



In this paper we studied the class of models in which the SM Higgs sector is extended by a 
single large complex scalar multiplet in such a way that the Higgs potential preserves an accidental 
global U(l) symmetry at the renormalizable level. The accidental U(l) symmetry is present when 
n = 2T + 1 > 5. Perturbative unitarity of weak-interaction scattering amplitudes involving the 
large multiplet excludes multiplets with n > 8. We choose the hypercharge to be Y — 2T so that 
the lightest member of the large multiplet can be electrically neutral. 

We showed that the models with n = 6, 7, or 8 are excluded by the incompatibility of the standard 
thermal freeze-out relic density with the dark matter direct detection cross section limit, assuming 
a standard thermal history of the universe. The model with n = 5 evades the direct detection 
constraint because its lightest state can decay via a Planck-suppressed dimension-5 operator during 
the first few days to few years after the big bang. 

This exclusion can be evaded if the model is modified in such a way as to break the global U(l) 
symmetry. One approach is to add one or more additional, smaller scalar multiplets in such a 
way that the accidental global symmetry is eliminated. This induces effective higher-dimensional 
operators involving the SM Higgs field and the large multiplet X that break the global U(l). Such 
operators typically also induce a vev for X and mixing between the states of X and the SM Higgs 
doublet. Models of this type involving a large multiplet with n — 7 and Y — 4 (whose vev preserves 
p = 1 at tree level [HI [27]) have recently been discussed in Ref. [5]. 
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A second approach is to arrange the model in such a way that the global U(l) is broken down to 
Zi. In this case, the real and imaginary components of the neutral member of X can be arranged 
to have different masses, so that the Z-mediated direct-detection scattering process x°' r N — > x°' l N 
is kinematically forbidden. The /i-mediated process x°' r N x°' r N has a much smaller cross 
section and remains experimentally viable. Such a model can be constructed for a large multiplet 
with n = 6 or n = 8 if its hypercharge is chosen as Y = 1. Then the U(l)-breaking operator 
($T T a (j>)(-X7T X) appears in the scalar potential, where $ and X denote the conjugate multiplets. 
We will address such models in Ref . [T2] . 
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Appendix A: Feynman rules 



We summarize the Feynman rules for the couplings of \ states to gauge and Higgs bosons. The 
following are for all particles and momenta incoming. We take Q > and denote the antiparticle 
of as X®* ■ 

Couplings to one or two Higgs bosons: 



X Q *x Q h = -iv 
X Q *x Q hh = -i 



A 2 - iUs ( Q - 



n - 1 



(Al) 



Couplings to two gauge bosons: 

,,2 



x Q *x Q w+w- = 



le 

■ o 

le 



(2Q + i)KA-Q 2 



9^u 



X (Q+2) *X Q W+W+ = — y/(Q + 2)(n - 2 - Q)(Q + l)(n - 1 - Q) g» v = X Q *X Q ^W~W~ 



x Q *x Q KK 
x Q *x Q z»z v 



2«e Q g^ v 
2ie 2 



X Q *X Q+1 W-A„ 



x Q *x Q+1 w-z v 

X Q *X Q Z,A lJ 



„2 2 
b W c W 

ie 
sw 



n- 1 



ie 2 / (Q + l)(n-l^Qj n 

? wV 2 (1 + 

ie 2 /(Q + l)(ra- l-Q) 



2Q)g^ = X (Q+1) *X Q W+A U 



SwC W 

2ie 2 Q 

SwC\V 



c^vQ 



2 

n-1 



[(1 + 2Q)c 2 w - n + 1] 9liV = x iQ+1) *X Q W+Z v 
9„u. (A2) 
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Couplings to one gauge boson: 
X Q *(p)x Q+1 (k)W- = 



ie_J(Q + l)(n-l-Q) 
s w 



( k -p)u = X 



- vW+i) 



(p) X Q (k)W+ 



X Q *(p)x Q (k)A ll = ieQ(k-p) L 



X Q *(p)x Q (k)Z, 



le 



SwCw 



n- 1 



(k~p) L 



(A3) 



Appendix B: Contributions to the oblique parameters from a scalar electroweak multiplet 

For a multiplet of hypercharge Y and size n = 2T + 1, the contribution to the S'-parameter is 
given by [H] 



S = 



4sf„c- 



a 



r'z Z (o) 



CwSw 



i%(o)-r; 7 (o) 



Y 
Gtt 



^Tflog(r 



(Bl) 



where m; and Xf denote the mass and third component of isospin of the complex scalar mass 
eigenstate Xi\j\ 

The contribution to the T-parameter can be similarly represented as |14j 



T = 



T W w(0) r zz (o) 



£m 2 log(m 2 ) [T(T+l)-(7; 3 ) 2 ] 



n— 2 



£(T - If )(T + 7? + l)/ 2 (m,,m, +1 ) 



where 



/a (nil, 7712) = I dx \xm\ + (1 — ar)m 2 ] log [xm 2 + (1 — ar)m 
Jo 



Finally, the contribution to the {/-parameter is 



U 



4s 



ii' 



[r' w (0) - c 2 w T' zz (0) - 2s w c w T' Zl {Q) - s 2 w T 



2 f 

77 V 



n— 2 n— 1 

£(T - T 3 )(T + 1? + 1) /iKnw+i) - - ]T(T?) 2 log(m 2 



i=0 



i=0 



where 



f\{mi,ra2) — I dxx(l — x) log [xm 2 + (1 — a^m 2 ;] 
j 



(B2) 



(B3) 



(B4) 



(B5) 



Note that we use the convention Q = T 3 + Y/2 and as such Y in the results of [14| must be replaced by Y/2. 
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We constrain these contributions to oblique parameters using a \ 2 variable including the corre- 
lations in the measured S, T, and U values, 



X 2 = °T P )(Oj - Of" )[<>*]£, (B6) 



where Oi is the ith observable and [a ]^ is the inverse of the matrix of uncertainties, 



[o*\ i j = kO i MO j p ij , (B7) 

where pij are the relative correlations (note pa = 1). For the three-observable case of interest, we 
can invert the matrix a 2 explicitly and write 



.2 



(! ~ PST - PtU ~ PSU + 2 PSTPTUPSu) 



(1 - p 2 TU ) (S - S cxp ) 2 (l-pl^iT-T^f 



(A^p) 2 (AT cxp ) 2 



(1 - p% T ) (U - C/exp) 2 . , AS- 5 cxp ) (T - T cxp ) 

* (PST - PTUPSU) 



(A[/ cxp ) 2 ^ u '^> AS exp AT cxp 

(B8) 



9 , x (T-Te.p)^-^) (5-5 exp )(C/-?7 exp ) 
(Ptu - PstPsu) r-^ ^ (Psu - PtuPst) 



A7^. X p A£/ cxp A^e^cj) AL/ exp 



Here S eX p) 2cxp, and f7 cxp are the experimental central values, AS cxp , AT cxp and At/ exp are their 
la experimental uncertainties, psr, Psu-, & n d Ptu & re their relative correlations, and S, T, and U 
are the contributions from the scalar multiplet computed using the formulas above. 



Appendix C: Contribution to h 



77 



The experimental observation of pp — > h — > 77 with a rate close to its SM value allows us to 
constrain the strength of the x°X *^ coupling in the Higgs pole annihilation region, m x o ~ rrih/2. 
The charged members of the X multiplet contribute to the loop-induced h — > 77 partial width [35] , 



T(h -> 77) = 



1024tt 3 



(CI) 



where i runs over charged particles of spin 1, 1/2, and 0, Q is the electric charge in units of e, N c . 
is the color multiplicity and the functions Fi(t) depend on the particle's spin, 



Fx = 2 + 3r + 3r(2 - r)/(r) 
F 1/2 = -2t[1 + (1 - r)/(r)] 
F = /3r[l-r/(r)]. 



Here r = 4mf/m^, and the function /(r) is given by 



(C2) 



f(r) = 



arcsm 



if r > 1 
if r < 1, 



(C3) 
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where we have defined r}± = 1 ± \f\--r. 

For the scalars, the coupling to the Higgs is parameterized by 



mf due to Higgs 



v 2 (X 2 /2-X 3 Tj/4) v 2 A Q /2 
AP+v 2 (X 2 /2~X 3 T l 3 /A) M 2 +v 2 A Q /2' 



where If is the third component of isospin of the scalar \i an d we define 

A Q = \2 + \\ 3 (n-l-2Q). 



(C4) 



(C5) 



Appendix D: Cross sections for co-annihilations 

The cross sections for annihilations of charged \ states into 77, Z7, W ± j, and // relevant for 
co-annihilations below the WW threshold are given by, 



OV re l{X Q *X Q -> 77) 

<?v rel { X Q *X Q -> ^7) 
^ rei ( X Q *X Q+1 -> ^ + 7) 



e 4 Q 2 



(n - 1) 



1 TO Z 

4m 2 , 



Q 



e 4 (Q + l)(n-l-Q)(2Q + l) s 
327rm x Q m x Q+i s 2 ^ 

= CTVre/ (x (Q+1) *x Q ^ 



Mr 



(m x Q + m x Q+i) 2 



Nc 

4-7T 



3/2 



f Q 



m x Q J ( 4w x Q - TO h) 2 ' 



where Aq was defined in Eq. ( C5 ) 
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